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Abstract 

The fractional Planck distribution is calculated by applying the Caputo fractional derivative with order p {p > 0) to 
the equation proposed by Planck in 1900. In addition, the integral representation of the Mittag-Leffler function is 
employed to obtain a new formula for the fractional BE distribution, which is then used to analyze the NASA COBE 
monopole data. Based on this analysis, an identity p ^ is found, where p is the dimensionless constant chemical 
potential that was introduced to the BE distribution by the NASA COBE collaboration. 
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1. Introduction 


There has been recent progress in physical sciences studies lEliiiBIS] that are based on fractional calculus 
iHiilBIIIllIl fisll . As part of these contributions, we investigate the NASA COBE monopole data llT3.fl5lfl^ 
by utilizing the Bose-Einstein (BE) distribution and a fractional calculus based distribution 117, 3 S. A well-known 
solution for Kompaneets equation 119, 3, which describes the photons distribution in the early Universe, is given 
by: 
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where x - hvIksT and Cb - 2hv^ jc^. yu is a dimensionless chemical potential. See Table[T]I) Kompaneets equation. 

On the other hand, the fractional calculus based photons distribution of the Universe is given by (See Table [T] II) 
Eractional calculus and III) Planck distribution in 1900 \ 2^ 2^ l^]): 
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where p is the fractional parameter. Ep(xP) is the Mittag-Leffler (ML) function defined as 
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The following digamma function i/r is used in ref. to analyze the NASA COBE data. 
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where i/r(z) = d(\nr(z))/dz. 

The COBE data is analyzed using Eqs. ([T]) and (|2ll, thereby yielding the following estimated values; 

T = l/kBj3 = 2.725 K, \ii\ < 7.58 x 10^^ 

\p-l\< 8.09x10^^ \p\xli < 5.47x10^"^ (1.00), 

1/7 - l|x/ 2 < 5.60 xlO^'^ (1.03), 

where, with Riemann’s ( function ^(3), 
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We would like to pay our attention to these similar values, p and /? - 1. The ratio between \p\ x I\ and \p - 1| x /2 is 
1.00 to 1.03. 

The same analysis of the COBE data is presented in this study, where we would like to adopt an integral represen¬ 
tation of the ML function: 
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A detailed derivation of the integral representation is supplied in ^ The magnitude of the integral representationi5(x, p)) 
contribution is estimated through concrete analysis of the COBE data in 0 Through such analysis of the COBE data, 
the following relation ensues; 


p - - \np. 

The concluding remarks and discussion are provided in ® 


(8) 


2. Integral representation of the Mittag-Leffler function 

The integral representation of the Mittag-Leffler (ML) function Eaiz) 1 11 1 ^ 12 , 13 ] is given by. 
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Table 1: Kompaneets equation, the fractional and the basic equations of the Planck distribution. 

I) Kompaneets Eq. describing the photons distributions in the early Universe 119. 201 
It is given by 

df/dt = C^x^^-§^x\{dfldxe + f + f) 

where Q = {kBTelmeC^){ne(Jelc) and Xe - hvIksTe- T^, T, n^, and cr^, are the electron temperature, 
radiation temperature, electron density, and Thomson’s cross-section, respectively. As a stationary 
solution, we derive the Bose -Einstein (BE) distribution with the chemical potential, 

fix) ^ l/(ce^ - 1) ^ - 1). _ 

II) Eractional calculus [71 

Eor the stationary solution in I) and a solution in III), if an inverse function Rix) - 1 ! fix) exists, 
the following equation is obtained, 
dRIdx - Rix) + 1. 

The Caputo derivative in fractional calculus is applied to the previous equation, i.e., 

QD^Rix) = Rix) + 1, 

The Caputo fractional derivative of the fix) function for m— 1,2,... is given by 
P‘’jix) = foix - t)'”-^-'/('”)(t)u't, (m - 1 < p < m), 

and 

lim oD‘’fix) = /"'Hx) = d"'fix)ldxf". 

p^m 

The following distribution is calculated using the ML function, 

fix)^l/Rix)^l/iEpixP)-l). _ 

III) Planck distribution in 1900 [21,22,23,241 

Planck utilized the following equation to describe the photons distribution U, 
dUIdp = -U - t/2. 

Adopting the aforementioned method in section II), the same expression is deduced. 


in the complex plane, z - x + iy i x - Re z, y - Im z), under the conditions. 


0 < a < 2, nal2 < 6 < min{7r, Trcr). 


( 11 ) 


As shown in Pig. [T] the contour 7 (e; b) (e > 0, 0 < 5 < tt) is comprised of the following three parts: (i) ray S si 
arg T = -b, |t| > e), (ii) circular arc Ca(0;e) ( -d < arg r < 6, |t| = e), (iii) ray 5a ( arg t = b, 1t| > e ). The left 
hand side of the contour ^(e; b), where the origin O is included, and is denoted by G^ fs', 6), while the other side is 
denoted by d)D 

After integrating along the contour •y(e; 6) while taking the e —» 0 limit into consideration, Eqs. (|9]l and (ITOl i 
respectively reduce to. 
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If 0 < a < 1, the explicit expression for Eaiz) is obtained through Eas. ([T2li and (ITsT l. If a > 1, using the ML 
function’s summation formula. 
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we can express Eaiz) through those with suffixes ajm < 1. Then, the ML function’s expression is given by the 
procedure in the case of 0 < a < 1. 
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Figure 1: The contour y(e; 6) composed of the ray S s, circular arc 0(0; s) and ray Ss 


In Fig. H we numerically examine the behavior of 6(p, x) in the ML function Ep{x'’) for Q < p < 2 and x > 0. 
The integral representation 6(p, x) for 0 < p defined in Eq.(|7|i approaches to zero in the limit of x —> oo. Function 
dip, x) at fixed p for 0 < p < 1 is negative and increases monotonously with x, and 6(p, x) for 1 < p < 2 is positive 
and decreases monotonously with x. For 0 < p < 2, the following limit can be derived through analytic calculations 

ill: 


|<5(p, x)|<|d(p, 0)1 = Ip-1 I/p. 


(15) 



p-1 

Figure 2: a) Contour maps of 6ip, jc) = constant for 0.5 < 



< 1. b) Contour maps of 6{p, x) = constant for 1 < p < 2. 
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3. Analysis of CORE data by Eqs. © and ® 


Expanding Eq. I©, providing that (p - 1) 1, as follows 
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= In p) - Ff^\x, 6(p, x)), 


we can analyze the NASA COBE data in Fig. [3] For comparison, the results on Ff^\x, In p) and F''^'^\x, 6(p, .x)) 
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Figure 3: Analysis of the COBE monopole data by Eqs. (T6), (T) and j4). a) Eq. {16). b) Eqs. 03- 0 and (T^. The magnitudes of error bars mean 
that 400 times real error bars. 


are shown separately. As is seen in Fig. [3] it can be said that contribution from is much bigger than that 

of Ff‘^\x, 6(p, x)). Then we can directly compare Eg. {p w ith e'/p and obtain Eq. (I8]i: p - -Inp. From this 
expression, p can be regarded as an inverse fugacity 126112711 . provided that NASA COBE data are used and the 
magnitude of 6{p, x) is very small. 

Analysis of NASA COBE data by Eq. (fThl) is presented in Table |2 


Table 2: Numerical proof of equivalences between Eqs. 0 and GD is shown at typical observed values of NASA COBE data. Our result of 
analysis by Eq. 1161 is as follows: T = 2.725008 ± 0.000026 K, |p - 1| = 8.09 X 10“^ and^^/NDF = 45.01/41. It should be noticed that we obtain 
the following values, T = 2.725016 ± 0.000008 K andx^/NDF = 45.10/42, provided that the first term is only utilized. 
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0.4.263569862x10-3 

2.920132007x10-3 

9.273229351x10-4 
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1.310610486x10-4 

4.892780012x10-^ 

1.016807344x10-'^ 


4. Concluding remarks and discussion 

Cl) We have shown that the ML function of Eq. (|6]l is decomposed into two functions jp and 6{p, x): The 
various properties of the integral representation 6{p, x) are investigated. In particular, through the analysis of the 
COBE monopole data, it is shown that the magnitude of 6{p, x) is much smaller than that of the Planck distribution. 
See, Fig. |3^). 

C2) In our previous paper |@1, we observed |p| = |p - 1| through the numerical analysis of the COBE monopole data 
by Eqs. ([T]i, (|2|i and Q. In the present paper, we have obtained Eq. (I8]l i.e., p = - Inp in the analytic form. 

C3) Combining the analytic relation of Eq. (I8]l with the numerical results of Table|2 the following relation is obtained, 

|p|=ln(|p-l| + l). (17) 
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D1) Using the present study, the interesting result given by Eq. (fTTI) is achieved. In the near future, the effect of the 
fractional parameter on other fields (BE condensation) should be investigated. Erom Eqs. (|2]i and we obtain the 
following formula. 


gyip) 
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It is well-known that this formula is given in Refs. Il26ll27ll28 l 


( 18 ) 
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